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Abstract 

We study double Higgs production in photon-photon collisions as a probe of the new dynamics 
of Higgs interactions in the framework of two Higgs Doublet Models. We analyze neutral Higgs 
bosons production and decay in the fusion processes, 77 — > SiSj, Si = h°,H°,A , and show that 
both h°h° and A A production can be enhanced by threshold effects in the region E^^ ~ e Zni}j -j- . 
Resonant effects due to the heavy Higgs, H°, can also play a role in the cross section enhancement 
when it is allowed to decay to two light CP-even h° or to two light CP-odd ^4° scalars. We have 
scanned the allowed parameter space of the Two Higgs Doublet Model and found a vast region of 
the parameter space where the cross section is two orders of magnitude above the Standard Model 
cross section. We further show that the Standard Model experimental analysis can be used to 
discover or to constraint the two Higgs doublet model parameter space. 

PACS numbers: 12.60.Pr, 14.80. Cp 
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I. INTRODUCTION 



The search for Higgs bosons is the prime task of CERN's Large Hadron Collider (LHC), 
with operation scheduled now for 2009. With the LHC guidance, the International e + e~ 
Linear Collider (ILC), which is currently being designed, will further improve our knowledge 
of the Higgs sector if that is how Nature decided to create mass. It was demonstrated in 
Ref. [l| that physics at the LHC and at the ILC will be complementary to each other in 
many respects. In many cases, the ILC can significantly improve the LHC measurements. If 
a Higgs boson is discovered, it will be crucial to determine its couplings with high accuracy, 
to understand the so-called mechanism of electroweak symmetry breaking ^]. The high 
resolution profile determination of a light Higgs boson (mass, couplings, self couplings, etc.) 
can be carried out at the ILC, where clear signals of Higgs events are expected with back- 
grounds that can be reduced to a manageable level. This is exactly the case of processes such 
as e + e~ — > 77 — > SiSj where Si = h°, H°, A . This fusion process can produce a Standard 
Model (SM) Higgs boson or one predicted by the various extensions of the SM, such as the 
Minimal Supersymmetric Standard Model (MSSM) or Two Higgs Doublet Models (2HDM). 

According to its Reference Design Report the ILC will run at an energy of yfs = 500 
GeV with a total luminosity of £ = 500/6 -1 within the first four years of operation and 
£ = \000fb~ 1 during the first phase of operation with y/s = 500 GeV. An e + e~ collider 
is uniquely capable of operation at a series of energies near the threshold of a new physics 
process. This is an extremely powerful tool for precision measurements of particle masses 
and unambiguous particle spin determination. Various ILC physics studies, indicate that 
a y/s = 500 GeV collider can have a great impact on understanding new physics at the 
TeV scale. An energy upgrade up to y/s ~ 1 TeV would probably open the doors to even 
greater discoveries. Another very unique feature of the ILC is that it can accommodate a 
77 collider with the photon beams generated by using the Compton backscattering of the 
initial electron and laser beams (J). In this case, the energy and luminosity of the photon 
beams would be of the same order of magnitude of the original electron beams. As the set 
of final states at a photon collider is much richer than that in the e + e~ mode, it would open 
a wider window to probe new physics beyond the SM. 

Since photons couple directly to all fundamental fields carrying electromagnetic charge, 
77 collisions provide a comprehensive means of exploring virtual aspect of the SM and its 
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a. 



extensions [5(. The production mechanism in hadron and e + e~ machines are often more 
complex and model-dependent. Thus, a 77 collider is much more sensitive to new physics 
even at higher mass scales (J. 

The primary mechanism of neutral Higgs boson production in 77 collisions is 77 — > 
flay , [loj], but in order to explore the triple and quartic Higgs couplings 
at future high energy colliders, it is necessary to study the Higgs boson pair production 



process. The trip 



linear colliders 



e Higgs couplings of the 2HDM have been extensively studied at e + e 



111 ] and shown to provide an opportunity to measure those couplings. At 



photon-photon co 



calculated in 
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liders, the cross section for neutral Higgs boson pair production has been 



13[ in the SM and found to be rather small. In the 2HDM, the process 



77 — > h°h° has been computed in the decoupling limit in Q, [jjj]. They found that the cross 
section can be substantially enhanced in the 2HDM and that the number of events expected 
at the Photon Collider will allow a determination or exclusion of some of the parameter 
space in the 2HDM potential. 

In the MSSM, various studies for Higgs pair production at a photon collider have 
been performed. The process 77 — > h°h° was studied in while reactions 77 — > 
h°H°, h°A°, H°H°, H°A° were determined in Il7|. Process 77 -> A°A° was calculated for 
the MSSM 18|, [l9| and shown to have a cross section of the order of 0.1 — 0.2 fb for a vast 



range of the photon-photon center of mass energy. 

In this paper, we present a complete calculation of pair production of all neutral Higgs 
bosons at the one loop level in the 2HDM. We study the Higgs self couplings effects on 
the 77 — > h°h° and 77 -> A A cross sections and briefly comment on the 77 -> h°A°, 
77 — > h°H°, 77 — > H°A° and 77 — > H°H° production modes. This exhausts all possible 
neutral scalar production processes in the 2HDM. A measurement of these processes can shed 
some light on the 2HDM triple Higgs couplings. However, even if the situation regarding 
a measurement of the vertex is not clear because no peak is detected, a vast region of the 
2HDM parameter space will be excluded. The scalars will be detected via similar final states 
because both the h° and the A , when not too heavy, decay predominantly into fermions. 
In this regard, the knowledge of their exact total cross section and angular distributions 
may be helpful in order to distinguish between CP-even and CP-odd scalars. Moreover, 
it is well-known that in the 2HDM, both the CP-even h° and the CP-odd pseudo-scalar 
A can be rather light 20[. In fact, the bounds on the h° and A masses originate from 
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the e + e~ — > h°Z and e + e~ — > h°A° production processes with the Higgs decaying to some 
combination of jets (mainly b jets) and r leptons. The production process e + e~ — > h°Z is 
proportional to sin 2 (a — (3) and this is the reason why LEP does not limit the mass of a 
light Higgs h° for sin(a — (3) = 0.1. For sin(a — (3) = 0.3 the bound is of the order of 80 
GeV [2l|. The pseudo-scalar mass is only limited by the results on e + e~ — > h°A°. However, 
if the sum of the masses is above the LEP energy limit, again no bound applies. 

A very interesting feature of 77 — > A A is that a light A can easily emerge in the Next- 
to Minimal Supersymmetric Standard Model (NMSSM) and therefore comparison between 
models will certainly prove useful. In addition, we also take into account in our calculation 
the perturbativity, unitarity as well as vacuum stability constraints on the various parameters 
in the Higgs potential. We will show that after imposing those constraints, cross sections 
are still large enough, in the hundred of fempto-barn (fb) region in some cases, to probe the 
2HDM scalar sector. We will also study some of these processes in the decoupling limit and 
in the fermiophobic limit of the so-called type-I 2HDM. 

The paper is organized as follows. In the next section, we review the 2HDM potential we 
will be using, give the analytical expressions for the triple and quartic Higgs couplings and 
list the theoretical constraints on the 2HDM scalar potential such as unitarity and vacuum 
stability. In Section IHIt we evaluate the double Higgs production cross section, 77 — > SiSj 
with Sij = H°, h°, A , in the general 2HDM paying special attention to 77 — ■> h°h° and 
77 — > A A . We then proceed to Section ITVl where we present our numerical results for the 
general 2HDM and for two limiting cases: the decoupling limit and the fermiophobic limit. 
In Section |V] we discuss the final states in the different 2HDM scenarios. Our findings are 
summarized in Section [VTl 



II. REVIEW OF THE TWO HIGGS DOUBLET MODEL 



A. The Two Higgs doublet model 



Two Higgs doublet models are some of the most well studied extensions of the Standard 
Model. Various motivations for adding a second Higgs doublet to the Standard Model have 



been advocated in the literature 



22 



231 ] . There are several types of 2HDM. While the 



coupling to gauge bosons is universal, there are many ways to couple the Higgs doublets to 
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matter fields. Assuming natural flavor conservation 24| there are four ways to couple the 



Higgs to the fermions 25] . The most popular models are the type-I and the type-II models, 
denoted by 2HDM-I and 2HDM-II, respectively. In 2HDM-I, the quarks and leptons couple 
only to one of the two Higgs doublet which is exactly what happens in the SM. In 2HDM-II, 
one of the 2HDM fields couples only to down-type fermions (down-type quarks and charged 
leptons) and the other one only couples to up-type fermions in order to avoid the problem 
of flavor-changing neutral currents (FCNC's) at tree-level. There are two additional models 
less discussed in the literature: model III in which one of the doublets couples to all quarks 
and the other couples to all leptons and the type IV model is instead built such that one 
doublet couples to up-type quarks and to leptons and the other couples to down-type quarks. 
There is a class of models sometimes called also type-Ill and denoted as 2HDM-III where 
FCNC are induced at tree-level 26( which can lead to fine-tuning issues. The discussion of 
the different models is not crucial as the production processes discussed here have a very 
mild dependence on the diagrams with fermion loops. It can however become relevant when 
discussing the different final states. When the Higgs decays predominantly to fermions, the 
relative size of the h° — > bb, h° — > cc and h° — > t + t~ branching ratios do depend on the 
model chosen regarding the Yukawa sector. Finally, we note that a 2HDM-I can lead to a 



fermiophobic Higgs boson h° 27| with suppressed couplings to the fermions (exactly zero at 
tree-level). In this case, the dominant decay mode for the lightest Higgs boson is h° — > 77 
or h° — > W + W~ , depending on its mass. No other version of the 2HDM possesses such a 
feature. 

The most general scalar potential, renormalizable, CP-conserving, invariant under 



SU(2) L ® U(l)y can be written as [22j: 



$ 2 ) = m?$!$i + m^l$ 2 + (m? 2 $t$2 + h.c) + ^A^i^) 2 + ^A 2 (^$ 2 ) 2 

+ A 3 ($ t 1 $ 1 )($ 2 $ 2 ) + A 4 ($ t 1 $ 2 )(^$ 1 ) + ^A 5 [($l$ 2 ) 2 + h.c] , (1) 

where $1 and $ 2 have weak hypercharge Y = 1 and vacuum expectation values (VEV's) v\ 
and V2, respectively, and Aj and m 2 2 are real- valued parameters. Note that this potential 
violates the Z 2 discrete symmetry $x — > $ l5 $ 2 — > — $ 2 softly by the dimension-two term 
m 2 2 (${$ 2 ), and has the same general structure as the scalar potential in the MSSM. From 
hermeticity of the potential, one concludes that m\ 2 must be real valued, so it can take both 
positive and negative values. 



After electroweak symmetry breaking, the and Z gauge bosons acquire their masses. 
Explicitly, three of the eight degrees of freedom in the two Higgs doublets correspond to the 
three Goldstone bosons (G , G°) and the remaining five become physical Higgs bosons: h°, 
H° (CP-even), A (CP-odd), and with masses m^o, ra#o, m^o, and m#±, respectively. 

The potential in Eq. (pQ) has ten independent parameters (including V\ and v^). The 
parameters m\ and m 2 are fixed by the minimization conditions. The combination v 2 = 
v\ + v\ is fixed as usual by the electroweak breaking scale through v 2 = (2\/2Gf)~ 1 - We 
are thus left with seven independent parameters; namely (Aj)i = i ; ... j5 , m 12 , and tan/3 = v 2 /vi. 
Equivalently, we can take instead 

m h o , m H o , m A o , m H ± , tan/3 , a and m 12 , (2) 

as the seven independent parameters. The angle (3 diagonalizes both the CP-odd and charged 
scalar mass matrices and a diagonalizes the CP-even mass matrix. One can easily relate the 
physical scalar masses and mixing angles from Eq. (CD) to the potential parameters, A.;, niyi 
and Vj. and invert them to obtain Aj in terms of the physical scalar masses, tan/3, a, and 
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29]. 



B. Theoretical and experimental constraints 

There are several important constraints on the 2HDM parameters imposed by experimen- 
tal data. In our analysis we take them all into account when the independent parameters 
are varied. 

First, the LEP direct search result in the lower bounds m^o > 114 GeV for a SM-like 
Higgs and m^o #o #± > 80-90 GeV for supersymmetric models in the case of the neutral 



scalars and for more general models in the case of the charged Higgs (see 21| for details). As 
stated in the introduction, the bound on the lightest CP-even Higgs heavily depends on the 
value of sin(a — f3). In a general 2HDM all bounds on the Higgs masses, with the exception 
of the charged Higgs, can be avoided with a suitable choice of the angles and m-^ 



Second, the extra contributions to the 5p parameter from the Higgs scalars 30J should 



not exceed the current limit from precision measurements [21] : \5p\ < 10~ 3 . Such an extra 
contribution to 5p vanishes in the limit m#± = m A o. To ensure that 5p is within the 
allowed range, we demand either a small splitting between m H ± and m A o or a combination 
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of parameters that produces the same effect. 

Third, the constraint from B — > Xg'-f branching ratio 



31 



321 ] gives a lower bound on the 



charged Higgs mass, m H ± > 295 GeV, in 2HDM-II. These bounds do not apply to model 
type-I and therefore are not taken into account in the fermiophobic scenario. Recent data 
from B — > Iv can also giv e a constraint on charged Higgs mass especially for large values of 



tan (3 in 2HDM-II 
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34]. 



Fourth, values of tan (3 smaller than « 1 are disallowed both by the constraints coming 
from Z — > bb and from B q B q mixing 31]. 

Finally, we should take into account the theoretical constraints. Let us start by noting 
that all 2HDM are protected against charge and CP-breaking 35J. We consider the per- 



turbativity constraints on the A, as well as the vacuum stability conditions [36( that assure 
that the potential is bounded from below. We require that all quartic couplings of the scalar 
potential Eq. [1] remain perturbative by imposing |Aj| < 8ir for all %. For the vacuum stability 
conditions we use those from 36J, which are given by: 



Ai > 



A 2 > 



VA1A2 + A 3 + min (0, A 4 - |A 5 |) > . 



(3) 



The above perturbative constraints are slightly less constraining than the full set of unitarity 



constraints 



37 



established using the high energy approximation as well as the equiva- 
lence theorem. It turns out, that requiring only perturbativity constraints on the X's could 
lead to scalar particles having a decay width which could exceed their mass. The problem 



is cured when we use the full set of perturbative unitarity conditions which are given by 



|o±|, |6±|, \c±\, \d±\, |e 12 |, \f±\, \g 12 \ < 8vr 



with 




c ± 
ei 

/+ 
9\ 



A 2 )±W(A 1 -A 2 ) 2 + -(2A 3 + A 4 )2 



i |(Ai + A 2 ) ± v /(A 1 -A 2 ) 2 + 4Al} , 

d ± = \ {(Ax + A 2 ) ± V /(A 1 -A 2 )2 + 4A|} , 
(A 3 + 2A 4 - 3A 5 ) , e 2 = (A 3 - A 5 ) 
(A 3 + 2A 4 + 3A 5 ) , /_ = (A 3 + A 5 ) 
92 = (A 3 + A 4 ) . 



(4) 

(5) 

(6) 

(7) 
(8) 
(9) 
(10) 
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These are very restrictive constraints on the allowed range of the parameter space. All values 
presented in the plots are consistent with all theoretical and experimental bounds described 
in this section. 

III. 77 SiSj, Sij = h°, H°, A IN THE 2HDM 
A. About the one-loop calculation 




FIG. 1: Generic charged Higgs and gauge bosons vertex like Feynman diagrams for neutral Higgs 
production 77 — ► SiSj in 2HDM. In the figures (f> = h° or H°. 

All processes 77 — > SiSj, Sij = h°,H°,A° with neutral Higgs in the final state are 
forbidden at tree-level and are mediated at one-loop level by vertex diagrams as well as 
by box diagrams. All those processes are sensitive to virtual gauge bosons, fermions and 
charged Higgs particles. We display in Fig. [1] and in Fig. [2] the generic Feynman diagrams 
with charged scalar particles exchange that contribute to 77 — > SiSj processes. Note that 
in Fig. [1] and in Fig. [2] we do not show the SM contribution with fermions and gauge bosons 
exchange. We have checked that they are always negligible when compared to diagrams 
with scalar exchange. The later comprise onedoop photon-photon fusion diagrams resulting 
in <fi = h° or H° intermediate states, followed by the decay h°, H° — > SiSj with Sij = h°,A° 



S 



(M (62) (63) (64) (6 B ) 




(6 6 ) (67) (6s) (69) (610) 




(611) (612) (613) 



FIG. 2: Generic charged Higgs and gauge bosons box like Feynman diagrams to neutral Higgs 
production 77 — ► SiSj in 2HDM. 

as shown in Fig. (TJ Ui_>4, vg and v±2- This kind of topology is sensitive to the triple Higgs 
couplings h°SiSj, H°SiSj, h°H + H~ and H°H + H~. For the 77 — > h°A° process, we have 
diagrams like i>i_»4, Vg and v i2 but with (p = A and also the contribution from the vertices 
with s-channel exchange of a Z boson v^ s , ^10 an d Vn. Note that for topologies like t>i_»4, v 9 
and W12 in Fig. [I] we have included the total width of the scalar particle <j) in the calculation 
of the corresponding amplitude. 

For the production mode 77 — > h°h°, the box contributions with virtual charged Higgs 
exchange is sensitive to the triple Higgs coupling h°H + H~ . On the contrary, in the case of 
the 77 — > A A process and again due to the CP nature of the pseudo-scalar Higgs boson 
A , it turns out that the box diagrams for 77 — > A A are rather sensitive to the A°H + G~ 
coupling which does not have neither a m 12 nor a tan/3 dependence. As one can see from 
Fig. [2] (diagrams (bi) and (611)), there are other topologies that contribute to 77 — * h°h° 
and 77 — > A A and which are sensitive to quartic couplings of the Higgs boson such as 
h°h°H + H~ and A°A°H + H~. 

As stated before, we are mainly concerned with the production modes 77 — > h°h° and 
77 — > A A . We will briefly comment on the h°A°, H°A°, h°H° and H°H° production pro- 
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cesses. The one-loop amplitudes were generated and calculated with the packages FeynArts 



391 ] and FormCalc 40j. The scalar integrals were evaluated with Loop Tools 4l|. The nu 



merical evaluations of the integration over 2 — > 2 phase space is done by the help of CUBA 



library 



421 ] . A cut of approximately 6° relative to the beam axis was set on the scattering 



angle in the forward and backward directions. 



B. Triple Higgs couplings 



The above processes are sensitive to triple and quartic Higgs couplings. Below, we list 
the relevant pure scalar couplings needed for our processes 77 — > h°h°, A A , h°A°. In the 
SM and in the general 2HDM these triple and quartic scalar couplings are given at tree-level 
by 



x sm 



h°h°h° 



\2HDM 
A h°h°h° 



\2HDM 
A H°H°H° 



\ 2HDM 
A H°h a h° 



\ 2HDM 
A H°H°h° 



\2HDM 
A A°A°h° 

\ 2HDM 
A A°A°H° 

\2HDM 
A AOG°hO 

\ 2HDM 
A H±HThP 

\2HDM 
A H±HTH0 



\2HDM 
A hOh°H-H+ 



\2HDM 
A A°A°H- H+ 



-3gm 2 h0 
2m w 

m w s 2 w 
-3<7 



\SM 

A h h°h°h 



-3g 2 m 2 



h° 



Am 2 w 



{cpc 3 a - spsl)s 2 pm 2 h0 + cl_ a c p+a m 2 l2 



m w s 2j3 
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2 m w s 2 



{cpcl - s p sl)s 2 pm 2 H0 + sl_ a s p+ 



a m l2 



2/i 



(2m 2 h0 + m 2 H0 )s 2a s 2 p + (3s 2a - s 20 )m 
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2/ii 



m 



h° 
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-g 
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g 
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\2HDM 
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2 m w 
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/ 3 1 3\ 2 1 2,2 
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m 



H± 



g 
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(11) 
(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
(19) 
(20) 



2m H ±s 2 pSp_ a + m h o(2c a +p + s 2a sp_ a )(2c a+ p — s^sp-a) 



C/3-a s 2a {cfi- a s w - 2s a+fi )m H0 + 



mTolc; 



2 2 



9 



2m w s 2/3 



2 1 _2 



m H0 (cf3_ a s 2f 3 - 2s{3 +a ) + m h0 {2cp +a - s 2 psp- a y 



(21) 
(22) 
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where g = e/ sin Q w is the SU{2) L gauge coupling constant. Here we use the short-hand 
notations sq and cq to denote, respectively, sinO and cos# where 9 stands for a or (3. All 
these triple Higgs couplings have a strong dependence on the physical masses (0 = 
hP,H°,H ± ,A°), on the mixing angles a and (3 and finally on the m 12 parameter which 
parameterizes the soft breaking of the Z 2 symmetry. 



ca 




100 200 300 400 100 200 300 400 

m 12 (GeV) m 12 (GeV) 



FIG. 3: The allowed regions for (1(77 — ► h°h°) in two Higgs doublet model. We have chosen 
m h o = 115 GeV, m A o = 270 GeV and m H ± = 350 GeV. On the left panel, m H o = 2m h o, 
i? 77 = 500 GeV, — 1 < sin a < 1 and 1 < tan/3 < 10. On the right panel, tan/3 = 1, sin a = —0.9 
and ELy = 800 GeV 



IV. NUMERICAL RESULTS 
A. The general 2HDM 

Before discussing our numerical results, it is worth pointing out that the following results 
are valid for all Yukawa type of couplings that do not generate FCNC at treedevel, as long as 
tan/3 remains small (tan/3 < 7 in the regions probed), as imposed by unitarity constraints. 
Moreover, as we will see later, the 2HDM contribution is dominated by scalar loops rather 
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than by fermion loops and the former are Yukawa model independent. Since data can easily 
accommodate light h° and A scalars [^(J in the 2HDM, we will concentrate hereafter on the 
h°h°, A°A° and h°A° production modes. In our numerical analysis, we will use: m t = 171 
GeV, m b = 4.7 GeV, m z = 91.187 GeV, m w = 80.45 GeV, the Weinberg angle s w is 
defined in the on-shell scheme as = 1 — m^/m|. For the fine structure constant we use 
a = 1/137.035989. 

We first note that we have reproduced the SM result for 77 — > H°H° and found perfect 
agreement with J12I . l^ . In the 2HDM case, our result agree with [3] while we have a full 
agreement with flji] if we take a = 1/128. The very detailed parton-level study [l^l conclude 
that for a 350 GeV center of mass energy photon collider and a Higgs mass of 120 GeV, an 
integrated 77 luminosity of 450 fb~ l would be needed to exclude a zero trilinear Higgs boson 
self-coupling at the 5a level, considering only the statistical uncertainty. If one assumes the 
luminosity based on the TESLA design report j^J we conclude that this is an attainable 
luminosity in approximately one year and certainly in less than two years. Therefore, we have 
decided to perform a comprehensive scan of the parameter space of the 2HDM looking for 
regions where the 2HDM dominate over the SM, that is, (J2Hdm{h — > h°h°) > gsm{h 
h°h°), together with perturbativity, unitarity and vacuum stability constraints on Aj. The 
results of this scan are shown in Fig. [3j In this scan we choose the charged Higgs mass to be 
350 GeV, in order to fulfill the b — > 57 constraint in all Yukawa type models. From the left 
scan, it comes out a correlation between tan (3 and myi. In order to have m\i large, unitarity 
constraints require tan f3 to be rather small. It is clear that in order to have a 2HDM cross 
section for 77 — > h°h° much larger than the corresponding SM one, a large m\i is needed 
together with a small value for tan/3 ~ 1, 1.5. With the set of parameters fixed, unitarity 
constraints forces tan (3 < 5 in the scanned region. The left scan shows that as tan (3 grows 
larger and larger, only smaller values of m 12 are allowed. In the right panel of Fig. [3] one 
can see that a quite large range of m# is allowed by perturbativity, unitarity and vacuum 
stability constrains. For large vn,\i > 260 GeV and for all values of m# < 500 GeV, the 
2HDM cross section for 77 — > h°h° can be 50 times larger than the corresponding SM cross 
section. From the scans presented, we conclude that for a large m 12 and even if tan f3 ~ 0(1) 
we have a significant slice of the parameter space where &2Hdm{i1 ~^ h°h°) can be much 
larger than the corresponding SM cross section while complying with all constraints both 
experimental and theoretical. 
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FIG. 4: SM and 2HDM total cross sections (7(77 -> h°h°) (left) and (7(77 -> A A ) (right) as a 
function of the two photon center of mass energy for unpolarized beams. Also shown is the total 
cross section in two other situations: both beams with right polarization and the two beams with 
opposite polarization. The parameters are chosen to be m h o, m H o, m^o, m\2 = 115,130,100,400 
GeV, sin a = 0.6, tan/5 = 1 and mjj± = 250 GeV. 

In Fig. H] (left) we show the polarized and unpolarized cross section for 77 — > h°h° both 
in the SM and in the 2HDM. The cross section is amplified by the threshold effect when 
E^ n m 2rriH± = 500 GeV, corresponding to the opening of the charged Higgs pair channel, 
77 — > H + H~ . Near this threshold region, the cross section of the 2HDM is more than two 
orders of magnitude larger than the SM one. Note that the s-channel vertex contribution 
is suppressed for large values of the center of mass energy due to the s-channel propagator. 
At high energies, box diagrams will dominates over the s-channel vertex. Let us take h°h° 
production as an example. The value of the couplings for this set of parameters is 

\2HDM ^ o v \SM y2HDM ^ r> r v xSM 

A h°H+H- ~ 6 x A h°h°h° ' A h°A°A° ~ x A h°h°h° 

\2HDM ^ on v \SM \ 2HDM ^ i 7 c v \SM 

A H a H+H- ~ zu x A h°h°h° 1 A H°A°A° ~ J-'-O X ^h°/i°/iO 

\2HDM ^ 7 X SM \2HDM ^ -, fi r \SM 

A h°h°h° ~ ' x A h°h°h,o 1 A H°h°h° ~ i0 -° x A h,0h°h° 

\2HDM ^ \SM \2HDM ^ n r xSM / 9 o\ 

A h°h°H+H~ ~ A h h°h°h 5 • A A0A0H+J?- ~ u -° x A h,0h0h°h° ■ 

Now the largest contribution to the vertex comes from the H°H + H~ and H°h°h° couplings 
squared which are approximately 400 and 272 times larger than (Afo^o) 2 - Note that in 
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FIG. 5: The total cross section of c(77 — ► hh) (left) and 17(77 ~~ > A°A°) (right) as a function 
of two photon center of mass energy in the 2HDM. With m h o,m H o,m A o,mi2 = 120,200, 120,300 
GeV, sin a = —0.86 and tan/3 = 1 for different values of m#±. 

the case of A A there are no boxes with virtual exchange because the A°H + H~ cou- 
pling does not exist and consequently the vertex contribution will dominate over the box 
contribution. 

In the right panel of Fig. H]the CP-odd Higgs boson pair production 77 — > A A is shown. 
The conclusions are very similar to the ones for the h°h° final state. The total cross sections 
is dominated by vertex contributions for low energy. This dominance is amplified by the 
threshold effect when i? 77 « 2m^±. For high center of mass energies, box contributions 
dominate because they have those t and u channel topologies which are enhanced for large 
center of mass energies. Finally, we show that when the initial photons are both right 

handed (++) or left handed ( ) polarized, the cross sections are enhanced by more than 

a factor two as compared to the unpolarized case. In the opposite case, when the initial 
photons have opposite polarizations, (-1 — ) or ( — h), the cross sections are now suppressed 
when compared to unpolarized case. 

The box contribution has either the quartic vertex h°h°H + H~ or twice the triple vertex 
h°H + H~ which means that the cross section will have (A/j0#+#-) 4 dependence. When 
h°H + H~ and/or h°h°H + H~ are large enough the boxes could dominate over the vertices. 
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By tuning sin a one can make h°H + H~ and h°h°H + H~ large enough to allow for the boxes 
to dominate. This is exactly what is illustrated in the following Fig. [51 where we present 
the total cross section for 77 — > h°h° (left) and 77 — > A A (right) as a function of the two 
photon center of mass energy for different values of the charged Higgs mass. Once more we 
can see the enhancement when i? 77 m 2rriH± for 77 — > h°h° and the cross section can reach 
0.2 pbarn. For this specific scenario, the triple and quartic couplings are 

\ 2HDM ^ 10 v \SM \2HDM ^ w \SM 

A h°H+H- ~ 16 X A h°h°h,o 1 A hPA0A° ~ » X A h°h°h° 

\2HDM r*^ A \s \SM \2HDM ^ r, .. \SM 

A H°H+H- ~ 4 x A h°h°hO j A .F/°A°A ~ z x A h°h°h° 

\ 2HDM ^ iSM \2HDM ^ n v \SM 

A h°h°h° ~ -J X A h°h°h° J A H°h°h° ~ "J X A h°h°h° 

A h°h.0H+H- ~ iz x A h°h°h°h° 5 A A0A0//+/f- ~ u -° x A h°h°h°h l /4 J 

In the case of the 77 — > /i /;, mode, the total cross section is now fully dominated by box 
contributions both at low and high energies. This is because the couplings h°H + H~ and 
h°h°H + H~ which contribute to the boxes are relatively large. In the case of 77 — > 74 y4 
mode, for low energy, the total cross section is dominated by vertex diagrams because 
the triple couplings h°H + H~ and /i°y4°y4 are large. For high energies, where vertex are 
suppressed, the total cross section is dominated by the box contributions. Hence, at high 
energies, both for the 77 — > h°h° and for the 77 — > A A modes, the total contribution is 
dominated by boxes. 

In Fig. [6l we illustrate the sensitivity to the charged Higgs boson contribution for different 
center of mass energy for the h h and the A A modes. In both cases, the cross sections are 
enhanced for light charged Higgs for the reasons explained above and are suppressed after 
crossing the threshold for 77 — > H + H~ production, i.e., rtin± > E 17 /2. 

In Fig. [7] we show the total cross section for 77 — > h°h° (left) and 77 — > A A (right) as 
a function the heavy Higgs mass for several values of the center of mass energy. Once the 
center of mass energy is close to m H , one can see in both plots the effect of the resonance of 
the heavy CP-even Higgs. The difference between them is only due to the vertex diagrams 
with an intermediate heavy Higgs that then decays to h°h° or A A . In both cases, the cross 
sections can reach 0.1 pb near the resonance w m#o. 

For the h°H° or H°H° final states, the generic set of Feynman diagrams that contribute 
to those processes is almost the same. The main enhancement factor for the cross section is 
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FIG. 6: The total cross section 17(77 — > h°h°) (left) and (7(77 — > j4°^4°) (right) as a func- 
tion of the charged Higgs mass for different center of mass energies i? 77 in the 2HDM. With 
m h o,m H o,m A o,mi2 = 120, 240, 270, 350 GeV, sin a = —0.9 and tan/3 = 1. 




FIG. 7: The total cross section 17(77 ~~ > h°hP) (left) and (7(77 — ► A A ) (right) as a function of 
the heavy Higgs mass m H o in the 2HDM. With m h o, m H ±, m A o, m 12 = 120, 250, 150, 200 GeV, 
sin a = 0.9 and tan/3 = 1.5 for different values of the two photon center of mass energy E„ n . 
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FIG. 8: The total cross section a(e + e~ -> 77 ->• h°h°) (left) and a(e + e~ -> 77 -> A°^4°) (right) 
as a function of the heavy Higgs mass mjjo in the 2HDM with the same parameters as in Fig. [7J 

the virtual charged Higgs bosons exchange, particularly relevant near the threshold region 
-E 77 = 2rriH±- The only difference between those final states and the h°h° and the A A 
ones, is the absence of the H° resonant effect, since it can not decay neither to h°H° nor to 
H°H°. The situation is the same as in the SM. If the CP-even H° has a mass of the same 
order as h° the cross section of h°H° or H°H° could be of the same order of magnitude and 
may reach 0.1 pb. If the CP-even Higgs is heavy, phase space suppression occurs and the 
cross sections for h°H° and H°H° production are smaller. 

In the case of h°A° or H°A° final state, a quite different situation occurs. Due to the 
presence of the CP odd scalar A in the final state, the Higgs boson <j) in the s-channel vertex 
(Fig. [T]) i>i->4, vg and V12 must also be CP-odd. Hence, the processes 77 — > h°A°, H°A° do not 
proceed through an intermediate CP-even Higgs. This implies that there is no closed loop 
of virtual exchange of charged Higgs, both in the vertex and box contributions, which again, 
is the main factor of enhancement of the cross section for the h°h° and the A A modes. In 
the boxes, we have rather a mixture of charged Higgs and charged Goldstones in the loop 
(gauge couplings). We have checked that the values of the cross section of 77 — > h°A° are 
much smaller than the h°h° and A A ones. We have performed a systematic scan over the 
2HDM parameters with m^o = m^o in the range 100 to 160 GeV, and found that the cross 
section of 77 — > h°A° does not exceed 0.04 pb. Note that the largest value for the cross 
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section, 0.04 pb, is attained near the top quark threshold region £" 77 pa 2mt] away from this 
region the cross section drops below 0.01 pb. Similarly, we found that large cross section 
for 77 — > h°A° prefer rather small values of tan (3 < 3, | sina| > 0.5, large m 12 and also low 
center of mass energies y/s < 600 GeV. 

Finally, in Fig. [8] we show the total cross section for e + e~ — ► 77 —>■ h°h° (left) and for 
e + e~ — >• 77 — > A°A° (right) as a function the heavy Higgs mass for two center of mass 
energies E^. The total cross section is evaluated by convoluting the photon-photon cross 
section with the photon-photon luminosity spectrum taken from the CompAZ library 44]. 
CompAZ is based on formulae for the Compton scattering and provides the photon energy 
spectrum for different beam energies and the average photon polarization for a given photon 
energy. First, let us remark that again this cross section is large enough to be measured in a 
significant region of the parameter space. We can still see the heavy Higgs resonance effects 
but somehow softened by the photon spectrum. 
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FIG. 9: The differential cross section for 77 — > h°h° and 77 — ► A ^ with equally polarized 
photons and oppositely polarized photons for y/s = 500 GeV (left) and y/s = 1 TeV (right). The 
2HDM parameters are the same as in Figj5] and mn± fixed at 250 GeV. 



Let us now turn to the differential cross section for 77 — > h°h° and 77 — > A A . In Fig. 
we illustrate the differential cross for a center of mass energy of y/s = 500 GeV and y/s = 1 

TeV and for identically polarized (++), ( ) and oppositely polarized (H — ), ( — h) initial 

photons. As one can see from the left panel for y/s = 500 GeV, for both the h°h° and 
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the A A modes, the angular distribution for (++) and ( ) polarized photons is almost 

flat while for the (H — ) and ( — h) modes it has a parabolic shape with a rather small cross 
section except for 77 — > h°h° where in the region of —0.5 < cos8 < 0.5 the differential cross 
section is of the order of 0.1 fb. For ^/s = 1 TeV, Fig. M (right panel), it is clear that the 
angular distribution are very similar to the ones on the left panel with the cross sections 
more than one order of magnitude smaller. The shape of these differential cross sections 
clearly tell us that we will not be able to distinguish the CP nature of the particle on the 
basis of this angular distribution. 



B. Fermiophobic limit 

In the SM, where just one doublet couples to all fermions, each scalar couples to the 
different fermions with the same coupling constant. In a general 2HDM it is also possible to 
couple just one doublet to all fermions by choosing an appropriate symmetry for both the 
fermions and the scalars. However, the difference between the SM and the 2HDM is that 
now the couplings are proportional to the rotation angles a and f3. For instance, the lightest 
CP-even Higgs couples to all fermions as cos a/ sin/3 Qj^f- By choosing a — ±|, the lightest 
Higgs decouples from all fermions. Such a scenario, with the appearance of the so-called 



"fermiophobic" Higgs boson, arise in a variety of models [271 ] . The heavy CP-even scalar 
will acquire larger couplings to the fermions than the corresponding SM couplings. All the 
remaining scalars are not affected by this choice as they do not couple proportionally to a. 

In the situation where the Higgs-fermion couplings are substantially suppressed, the full 
decay width of the Higgs boson is shared mostly between the WW, ZZ and 77 decay modes. 
In this limit, for masses m h o < 100 GeV, the Higgs boson dominantly decays to photon pairs. 
Experimental searches for fermiophobic Higgs bosons at the LEP collider and the Tevatron 
collider have yielded negative results. Mass limits have been set in a benchmark model that 
assumes that the couplings hWW and hZZ have the same strength as in the SM and that 
all fermion branching ratios are exactly zero. Combination of the results obtained by the 



LEP collaborations 



45|, 



46, 



47J, 



481 ] using the process e + e — > hZ, h — > 77 yielded the lower 
bound rrih > 109.7 GeV at 95% C.L. 49J. In Run I, Tevatron has set lower limits on rrih 



by the DO and CDF collaborations which are respectively 78.5 GeV [50[ and 82 GeV [51], 



using the processes qq' — > V* — > hW, hZ, h — > 77, with the dominant contribution coming 



19 



tan(3=2, m 12 = 200 GeV 
tanp=3, m 12 = 100 GeV 
' tanp=5,m 12 = 40 GeV 
s tanp=7, m 12 = GeV 



Fermiophobic 




400 600 800 1000 1200 1400 
E^GeV) 



10 u 



10" 



10 



10" 



10" 



10" 



-2 



I m 12 = 300 GeV 



Fermiophobic 



= 200 GeV\ 




m 12 = GeV 



400 600 800 1000 1200 1400 



E^GeV) 



FIG. 10: The total cross section (1(77 — ► h°h°) as a function of the center of mass energy 
in the fermiophobic limit (sin a = 1). In the left panel with different values of tan/3 and mi 2 for 
mf l o,m H o,m j ^o,mH± = 100, 200, 100, 150 GeV. In the right panel, m/jO = 120 GeV, m#o = 
m H ± = rrijifi = 300 GeV and tan/3 = 1. 



from the W boson. Recently the CDF 52J and the DO 53| collaborations have improved 
their bounds which are now of the order of 100 GeV, close to the bound obtained by the 
LEP collaborations. It should be noted that all experimental mass bounds assume, in the 
fermiophobic limit, tan/3 ~ in the tree-level couplings to the gauge bosons. 

It is clear that in the fermiophobic limit the coupling H°h°h° (Eq. (1141) ) is directly 
proportional to m 12 , while the other couplings depend both on m 12 , tan/? as well as on 
m,H±- We show in Fig. [TOl the total cross section for 77 — > h°h° as a function of the two 
photon center of mass energy in the fermiophobic limit. The parameters m\i and tan/3 
are varied within the allowed range. It is clear that the cross section is enhanced for large 
m 12 and can reach 0.1 pb for large m 12 = 200 GeV (left) and m 12 = 300 GeV (right). 
The observed kinks are the top threshold at _E 77 = 2m t and the charged Higgs threshold 



E 



77 



2m 
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C. Decoupling limit 



A study of 2HDM in the decoupling limit reveals the case where all scalar masses except 
one formally become large and the effective theory is just the SM with one doublet - m h o << 
m$ where m<j, = m H o H ± (see 28j] for an overview). In this case, the CP-even h° is the 



lightest scalar particle while the other Higgs particles H°, A and H ± are extremely heavy. 
In 2HDM, the decoupling limit can be achieved by taking the limit a — > (3 — n/2. This 
means that the coupling of the h° to the gauge bosons, fermions and light Higgs, h° are the 
same as for the Standard Model h SM Higgs. Also, in the decoupling limit, the triple Higgs 
coupling X^o h o H o vanishes at tree-level, so that the heavy Higgs cannot contribute to the 
process 77 — > h°h° and the result is independent of the mass m H a. In the decoupling limit, 
the tree-level trilinear Higgs couplings take the form 

\ 2HDM ^ ySM / 9 c\ 
A h°h°hO ~ A h°hOhO, {^) 

X 2HDM ^ n 



K ° H+H 2^ 



w 



m 2 h0 + 2m 2 H ± + m\ 2 



\ 2HDM _ 9 \ 2HDM / ofi \ 

A h°h°H+H- ~ 77~ — A h°H+H-- l /D J 

It is clear that these couplings are independent of tan j3 as well. As one can see from the 
analytical expression of ^oh+h- an< ^ A hOh°H+H-i ^ ne charged Higgs and m\ 2 add construc- 
tively for m\ 2 > and destructively for m\ 2 < 0. At tree level, the enhancement of the 
cross section essentially depends on the size of the h°H + H~ and H°H~H + couplings. By 
taking these couplings as large as possible under the referred theoretical and experimental 
constraints, we obtain the best possible enhancement for the cross section 17(77 ~~ * h°h°) in 
the 2HDM for each m$ (m$ = m#o = m^o = mjj±). However, it is well known that those 
couplings of the CP-even h° which mimic the SM couplings get significant radiative correc- 
tions in the decoupling limit to which we refer to as non-decoupling effects. Several studies 
have been carried out looking for non-decoupling effects in Higgs boson decays and Higgs 
self-interactions. Large loop effects in h° —>■ 77, h° — > and h° — > bb have been pointed 



out for the 2HDM 
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55J and may provide indirect information on the Higgs masses and 
the involved triple Higgs couplings such as Xhon+H- for h° — > 77 and \h°H H°, -^0,40,40 an d 
Xh°h°h° for h° — > bb. 

The non-decoupling contributions to the triple Higgs self-coupling XhOh°h° have been in- 
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vestigated in the 2HDM in Ref . 56( , using the Feynman diagrammatic method. It has been 
demonstrated that the one- loop leading contributions originated from the heavy Hi ggs boson 
loops and the top quark loops to the effective h°h°h° coupling can be written as [56| 



\eff 
A 



3 K° 



711 



HO 



M 2 



711 



AO 



M 5 



h°h°h° 



127r 2 m? f 2 
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M 2 



67T 2 m 2 t0 f 2 
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H±. 



N c m\ 
2>TT 2 rn\ v 2 
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11 1 11 



mipV 



m% v 




^ ,(27) 



where M 2 = m 2 2 /(sin/3 cos/3), m $ and pi represent the mass of H°, A or H ± and the 
momenta of external Higgs lines, respectively, N c denotes the number of colors, and m t is 
the mass of top quark. We note that in Eq. ( 1271) m^o is the renormalized physical mass of the 
lightest CP-even Higgs boson h°. In the calculation of the 77 — > h°h° cross section in the 
decoupling limit, we replace the X^d h o h o coupling by its effective coupling given in Eq. (T2 



which corresponds, in this limit, to an effective 2-loop 2HDM contribution (see Ref. 56j for 
a detailed discussion). 
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FIG. 11: Cross sections for hPhP production in the decoupling limit with unpolarized photons. On 
the right we show the loop contributions to the total cross section as a function of m$ and for two 
values of mi2, and 200 GeV. On the left panel the cross section as a function _E 77 is shown for 
different values of m$ and m\i- The light Higgs mass is m^o = 120 GeV. 

As stated in the introduction, this process was studied in detail in Refs. Q, [3| in the 
decoupling limit. In fact, there is a one to one correspondence between our potential and 
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the potential used in [jjj] which relate A5 parameter of Ref [ijj] to our mi 2 by: 



m 



12 



-X5V1V2 (28) 



From this relation, one can see that A5 < (resp A5 > 0) correspond to our m\ 2 > 
(resp m\ 2 < 0). In [ljj], the non-decoupling effects have their origin in the h°H + H~ and 
h°h°H + H- vertex and are present already at tree level. In their notation, when A5 = 
(m.12 = in our notation) and at the same time the charged Higgs mass and the lightest 
Higgs mass are of the same order we get Xh°H+H- ~ ^h°h°h°- Hence, in this limit the 2HDM 
cross section is very similar to its SM counterpart. As m\ 2 falls to negative values, the cross 
section drops due to a cancelation between the m\ 2 and the m 2 H ± contributions until the 
m\ 2 term starts to dominate and the cross section increase again. For positive m\ 2 the cross 
section is enhanced relative to the SM one and the non-decoupling effects can be quite large 
as shown in [l4| for A5 < 0. Therefore, at the 1-loop level m\ 2 and the charged Higgs mass are 



the parameters that regulate the non-decoupling effects. In [15J], the authors have studied 
only the case of m 12 = and looked for non-decoupling effects in higher order corrections 
with origin in the vertex h°h°h° as described earlier. In this case the non-decoupling effects 
appear mainly for l arg e values of the masses. In this section we will combine both the effects 
of ref [ijj] and ref 15] and present the results for the case of the unpolarized photon cross 
section. We will show that even in this case where the cross sections would be severely 
reduced, there are still regions where the 2HDM CP even Higgs h° could be disentangled 
from the SM h SM . 

In the left panel of Fig. [11] we show the cross section for 77 — > h°h° as a function of i? 77 
for m$ = 300, 400 and 500 GeV with m\2 = together with the case where vii\2 = 200 
GeV (note that from now on we will be considering m\ 2 > 0) and m$ = 300 GeV. In this 
left panel, the coupling h°h°h° is taken at the tree-level without the higher order correction 
of Eq. (l2Tj) . The non-decoupling effects due to the charged Higgs mass can be seen in the 
vicinity of -E 77 = 2 m$ where the cross section is of the order of 0.6 to 0.7 fbarn. This effect 
is enhanced for higher values of myi- As shown in the plot, for a charged Higgs mass of 300 
GeV and m\2 = 200 GeV the cross section can reach 5 fbarn and this value grows with large 
vrt\2. In the right panel of Fig. [11] we display the cross section of 77 — > h°h° as a function 
m<j,. Here, besides the SM value we plot four different scenarios. The one-loop case with 
rrt\2 = and mu = 200 GeV and the two-loop case with the higher order corrections given 
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in Eg. (1271) for the same two values of mi 2 . One can see that the cross section enhancement 
due to the large corrections in \% h0 take place only for large if m 12 = 0. As m 12 grows 
the cross section grows as described in the left panel, but on top of that we get an extra 
enhancement due to the higher order corrections. Largest values of the cross section, that 
can reach 10 fbarn, are attained for the low mass region in m$. It is also clear from the right 
panel that without the higher order corrections of Eq. (l2Tj) . for large m<j>, the 2HDM cross 
section is similar to SM one. On the other hand, once those higher order corrections are 
included, non decoupling effects arise - the cross section then grows relative to the SM one 
and can be almost two orders of magnitude above the corresponding SM process depending 
on the value of m$. We finally note that the cut on m$ at 610 GeV for m 12 = and on m$ 
at 550 GeV for m 12 = 200 GeV is due to unitarity constraints. 



V. HIGGS SIGNATURES 



We are considering a light CP-even Higgs, that is, with a mass of 120 GeV or less. 
Assuming that all decay channels with some other Higgs boson in the final state are un- 
accessible, this particle decays predominantly to bb in this mass region. The exception is 
in the fermiophobic Higgs scenario where it decays to two photons although for a mass of 
120 GeV one has already to consider the decay to two W bosons even if one of the W is 
off-shell and strongly virtual. The rate at which it decays to each final state depends on the 



remaining parameters of the 2HDM (see 
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for details). The two subleading decays 



that compete with h° — > bb are h° — > cc and h° — > t + t~. In model type I the branching 
fractions to fermions are the SM ones because the coupling dependence cancels. In model 
type II, the ratio T(h° — > bb)/T(h° — > t + t~) is the SM one. On the other hand it is easy to 
check that for tan (3 > 1 the decay h° — > bb is again the dominant one provided we are in a 
region with moderate values of tana. For the remaining Yukawa models the situation does 
not change dramatically. However, if we take as an example the case where $ 2 couples to 
the quarks and $1 couples to the leptons, we obtain the ratio 

rfo -> bb) mf 1 , , 

T(h°~^T+T-) m 2 T tan 2 a tan 2 /5 1 ' 

in the limit rrih 3> m q . Even for tana = tan/? « 3 the ratio becomes almost 100 times 
smaller than the corresponding SM ratio. Therefore, a detailed study for each model will 
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have to take into account the exact branching fractions for each Yukawa version of the 



The dominant background to double Higgs production is 77 — > W + W~ and non-resonant 
four jet production. The first one can be reduced by imposing a cut on the invariant mass of 
each pair of b-jets, M(qq), forcing it to be close to Higgs mass. An efficient b-tagging would 
further reduce the background by asking that at least three jets be identified as originating 
from b quarks. A cut on the polar angle would eliminate the non-resonant 4-jet background. 
Together they would reduce the backgrounds to a level well below the signal. A more detailed 
analysis can be found in [k]]. In the fermiophobic case, the analysis is greatly simplified by 
the smallness of the four photon production cross section. We just need to avoid very soft 
photons which can be done with a sensible cut on the photon's transverse momentum. 

The process 77 — > h°H° gives rise to very similar signatures. There are in principle two 
drawbacks: first the phase space is reduced because there is a heavier Higgs in the final state; 
second we can not reduce the background by asking the two invariant masses from each pair 
of b-jets to have a similar magnitude. On the other hand, if the channel H° — > h°h° is open, 
it can be dominant. This would lead to a very interesting signal of a six b-jet final state. As 
for 77 — > h°(H°)A° we concluded that the cross section is much smaller due to the absence 
of virtual charged Higgs in the loop. 

A pseudo-scalar Higgs decays again mainly as A — > bb and the subleading competing 
decays are again A — > cc and A — > t + t _ . The situation is similar to the CP-even Higgs 
case - if tan (3 « 1 the bb channel is always the dominant and well above the others. Once 
more if we consider the model where one doublet couples to the quarks and the other couples 
to the leptons, we obtain the ratio 



in the limit rrih 3> m q . In this case, for tan/3 pa 3 the ratio becomes 10 times smaller than 
the corresponding SM ratio. Note that a similar situation can occur in the ratio between 
decays to down and to up quarks. Obviously, this is true as long as the other channels 
involving other Higgs are closed. The decays to either h°Z or W + H~ become dominant as 
soon as they are kinematically allowed. The last two cases could again lead to interesting 
final states which are easy to detect. As all the cases discussed refer to light Higgs, even 
if the Higgs to Higgs channel is open the Higgs decay width will in most situation be well 



2HDM. 



T(A° -> bb) 




(30) 



r{A° -> t+t 



m 2 tan 2 /? 
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below the GeV or a few GeV at most. 

A final word about the behavior of the cross section with the scattering angle. We have 
shown that if the A and h° masses are of the same order, and because in most models 
the possible final states are very similar, it will be very hard to distinguish a CP-even from 
a CP-odd state. In fact, even if one changes the polarization of the initial photons, the 
differential cross section does not distinguish clearly between the two cases except in regions 
where either the cross sections are too small to be measured or the angle is too small to be 
probed. 

VI. CONCLUSIONS 

We have calculated the total cross section for 77 — > S{Sj, Si = h°,A°, H°, in the frame- 
work of the 2HDM taking into account perturbativity, unitarity as well as vacuum stability 
constraints on the scalar potential parameters A's. All available experimental constraints 
were also taken into account. For the numerical study, we mainly focused on the h°h° and 
A A production modes. We have studied those processes in the general 2HDM, and in 
two other limiting scenarios: the fermiophobic limit, a scenario where the lightest CP-even 
Higgs decouples from the fermions and the decoupling scenario where this same Higgs re- 
sembles the SM Higgs boson. We have shown that, for both production modes, the most 
important contribution to (7(77 — > h°h°) and to 17(77 ~~ * A°A°) comes from the charged 
Higgs H ± diagrams and also from the diagrams with a resonant heavy Higgs that can decay 
as H° — > h°h° or H° — > A A . Since the cross section is dominated by the charged Higgs 
virtual exchange, it does not depend on the Yukawa structure of the model. 

We have shown that the cross section for 77 — > h°h° can be more than 100 times larger 
than the corresponding SM one in vast regions of the parameters space. The parameter space 
will easily be probed for the largest allowed values of tan/3, m\ 2 and | sino;|. A light charged 
Higgs, that is, below the collider center of mass energy, is preferred. A variable energy 
collider would be a good option to detect the heavy Higgs resonance. We have argued that 
knowledge of the charged Higgs effects may be crucial to understand the nature of the Higgs 
bosons if they are eventually found in future experiments at LHC and/or ILC. In the case 
of 77 — > h°h°, we have illustrated that even with a charged Higgs mass of the order of 
300 GeV, in agreement with b — > 57, the cross section can have a substantial enhancement 
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(clearly at least one order of magnitude above the SM). In case of 2HDM type I, where a 
light charged Higgs is not ruled out from b — > 57 one can have an even larger enhancement 
which could be 3 order of magnitude above the SM results. 

The analysis in Q shows that the SM H.ggs triple coupling con.d be probed at a linear 
collider. As described before, their analysis is mainly based on an invariant mass cut, on 
the identification of at least 3 jets as originating from b-quarks and on a the polar angle 
cut I costal < 0.9. We have showed that the inclusion of the new 2HDM diagrams do not 
change the angular distribution as the Higgs angular distribution remains almost flat, so 
that the same cut could be applied. Moreover, in the two Yukawa versions of the model, 
BR[h — > bb) is at least the SM one if not larger. Because the invariant mass cut is the 
same, the analysis can be applied directly to the 2HDM case. Therefore, when a complete 
experimental analysis is completed for the SM, it is ready to be used to constraint the 2HDM 
parameter space. This is one of the major advantages of this study. 

Although other regions give rise to higher cross sections, the very interesting case of the 
2HDM decoupling limit can also be probed at the photon collider. The importance of the 
sign of m\ 2 was studied in a more general context. Clearly, positive m\ 2 O 11 our notation) can 
lead to large non-decoupling effects. Non-decoupling effects can also appear due to higher 
order correction to the triple h°h°h° vertex. We have shown that non-decoupling effects will 
be more easily seen in the low and in the high regions as in the intermediate mass region 
the cross section is closer to the SM values. 

In the fermiophobic limit this process is complementary to the LEP production process 
as it grows with 77112- Most importantly, it can also probe a part of the parameter space 
that cannot be accessed at hadron colliders. When m 12 ~ the cross section vanishes at 
hadron colliders 59j. On the contrary, we have shown that in photon-photon collisions the 
cross section can reach a few fbarn for m 12 rs and tan f3 ~ 5. The region of low m 12 in the 
fermiophobic scenario will most probably not be excluded until we have access to a photon 
collider. 

Regarding the CP-odd Higgs, we have shown that, for the energies considered, only a 
light A will be probed at a photon-photon collider. For m^o > 150 GeV the cross section is 
virtually zero. We have also covered the mixed CP-even modes h°H° and H°H° final state 
and checked that they can have comparable cross sections to the h°h° mode if H° is not 
too heavy. For the other modes, h°A° and H°A° we concluded that their cross section is at 
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least one order of magnitude smaller than the h°h° and the A A one. 

As for the final states we have shown that in the two most popular models, 2HDM-I and 
2HDM-II, the bb final state is the preferred channel as its branching fraction is always at 
least the SM one. In this case and for not too heavy h° and A a 46 final state can be 
searched for. We have showed the such a study for the SM can also be used for the 2HDM. 
This study can then be complemented with the final states 262r and At. In other models 
and when other channels are kinematically available, a detailed study has to be performed 
taking into account all 2HDM parameters. 

A last word: there is a well known complementarity between the LHC and the ILC In 
this regard, if some of the 2HDM parameters and/or couplings are measured at the LHC 
and/or ILC experiments such as Higgs masses and the magnitude of the vertices h°h°h° and 
H°h°h°, this information could then be used at 77 experiments in order to extract missing 
parameters and couplings like h°H + H~ and H°H + H~ . 
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